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^ . Abstract: 

^ ■ A beautiful and intriguing relationship has recently been proposed to express the critical screening lengths associated with the 
, apparition of new bound states for the two-dimensional statically screened Coulomb potential. Semiclassical quantum theory 
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show that this relationship is unfortunately not strictly exact. These results are confirmed by a variational calculation, which 
provides upper bounds for the critical screening lengths in the case I > 1. 



O ' I. INTRODUCTION 

The two-dimensional statically screened Coulomb potential Vg (r) plays a central role in the physics of semiconductor 
heterostructures near the band gap, especially with respect to their linear and nonlinear optical properties when the 
influence of excitons is properly taken into account Although already considered three decades ago it 

has nevertheless received a lot of recent attention and been analyzed using different approaches such as the WKB 
approximation L perturbation theory Q , variational calculation or numerical resolution based on the variable- 
phase method [|6f 7| . 

Using a variable-phase approach Q to count the bound states allowed by Vs{r) — only a finite number of them 
fH , exist — and making a connection with the Levinson's theorem in two dimensions, Portnoi and Galbraith recently 
Q ■ proposed a beautiful and intriguing relationship between the critical screening lengths (scaled to the exciton Bohr 
O radius) and the number of bound states for a given value of the angular momentum (Z > 0): 



(the original notation used v ~ + 1, whereas we focus on the number A^ of nodes of the wavefunctions); they 
admitted that its analytical derivation had still to be found They then proceeded to derive a formula giving the 
number of bound states as a function of rg , which markedly differs from a WKB estimate |^ and from the Bargmann 
bound condition reformulated for the two-dimensional case 0. 

In this work, we show that semiclassical quantum theory simply demonstrates that Eq. (|l]) is unfortunately not 
. strictly exact for I = states, and that a previous estimate using the WKB approximation is incorrect. Using a 
variational calculation, we also give upper bounds for {rs)^ in the case I > 1, which are nearly given by Eq. (|l|). 

■ ; II. SEMICLASSICAL QUANTUM THEORY AND THE TWO-DIMENSIONAL SCREENED COULOMB 
^ '. POTENTIAL 

O . . 

O ' The two-dimensional screened Coulomb potential Vs has a very simple form when expressed in wavevector space 



VM = ~^, (<ls^-) (2) 

q + qs \ TsJ 

where we have adopted the same units as in Ref. |^ (the length scale is the effective Bohr radius oq and the energy 
scale is the three-dimensional exciton energy R). In real space, it reads [^J^J^ 

^.(r)--^{l-|g.r[Ho(g.r)-A^o(gsr)]} , (3) 

where Hq and Nq are the Struve and Neumann functions, respectively. The term between braces in Eq. (^ goes 
down to like /r^ as r goes to infinity |^ , which is a much reduced decrease with respect to the exponential decay 
occurring in three dimensions. 
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Let us now turn to the semiclassical quantum theory approximation for a potential V{r) . Using conformal mapping, 
Yi and collaborators have shown that the two-dimensional formulation of energy quantization may be written as 

^ ' v/S-14ff(r) dr =(^N+^^ TT, (iV > 0) (4) 

where Veg{r) = V{r) + P /r'^ , while ri and r2 are the classical turning points. They have proved that Eq. ^ gives 
exact results in several cases, and excellent ones for the energy levels of impurity states in an arbitrary external 
magnetic field. Since we are looking for the critical values of at which new bound states appear at zero energy, we 
can set i? = in Eq. so that, for s eigenstates of Vg, 

y^{l-|9.r[Hofe r)-No{qs r)]} dr = (^N +^ t: . (5) 

In the following we call {rg)^ the value obtained by using the semiclassical approximation in contrast to the exact 
one. Eq. (||) gives 




I-^ ^ [^^{u^)-No{u^)\du. (7) 

As is usual with semiclassical expressions, we expect Eq. ^ to be accurate when N ^ \. If Eq. were true, we 
would find 2 J = tt, so as to get the correct leading dependence of {rs)c in N"^ ■ ^ has been numerically evaluated 
by first calculating the integral between and 10, and then bracketing the integral from 10 to +oo by using the 
asymptotic expansion of Hq — A^o- The same procedure, repeated for an intermediate bound of 20, gave the same 
result, namely 

21 w 3.14057414. (8) 

It definitely differs from tt, even though the relative error is only about 3.24 x 10"*. Consequently, for ^ = states 
and large iV's, (r^)^ should be different from Eq. (|l|). 

Portnoi and Galbraith ||^ compared Eq. (Q) with the WKB prediction of Reyes and del Castillo-Mussot for Z = 
states ||], and observed a disagreement by a factor 2.5. By contrast, Eq. (^) satisfactorily removes the discrepancy. 
Why then do our semiclassical approximation and the WKB result differ by such a large amount ? The explanation 
probably comes from the tricky use of WKB wavefunctions and asymptotic expansions in Eqs. (7)-(12) of Ref. 
In particular, their Eq. (12) contains w^"^/^ ^1 — w dw, which diverges. Admittedly, this integral can formally be 
written as B{3/2, —1/6), whose true value is -6.72 instead of the given -2.81. 



III. UPPER BOUNDS FOR THE CRITICAL SCREENING LENGTHS 

At this point, we have merely shown that Eq. cannot be exact for I = states, since it does not verify the large 
N limit. Strictly speaking, we cannot use this argument to reject the hypothesis that (rg)^ =1, 3, 6, 10, ... for / = 
and A' = 1, 2, 3, 4 ... or more generally, for low-energy states. However, because of the near-coincidence of 21 with 
TT, one may have second thoughts about the accuracy of {rs)c- K is very long range because of its 1/r^ decrease at 
infinity. Others calculations |y, not detailed here, seem to indicate that (rs)^(A — 0, 1 > 2) > I {2 1 + 1) . 

For this reason, we have also relied on a variational method to provide further information. Rescaling the lengths 
to X = r/vs, the Schrodinger equation verified by the radial wavefunctions V'o(a^) of zero energy reads 

<(a;) + ^ ^i(x) + {l - I X [Ho(x) - Noix)]] - M^) = 0. (9) 

We can then follow Hulthen and Laurikainen's method jl^ — first developed to determine critical screening lengths 
for the Yukawa potential in three dimensions — to calculate upper bounds for (rs)^ when I > 1 ||ll| . The principle of 
their method is to find extrema in A of J — A A, now defined by 
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,10) 

dxcP^{x)^(l-^x[IIo{x)-No{x)i^ (11) 



in which is a trial function expanded on a set {(pi^ ■ • • , <^m}- 

There are well-known difficulties with this approach: (i) the trial functions must be chosen so as to converge 
towards the true eigenfunction of Vg (ii) convergence rapidly deteriorates with N (iii) the number of significant digits 
with which the matrix elements of J and N must be computed greatly increases with M. None the less, preliminary 
numerical computations ||^ with admittedly nonoptimal {fj} and M = 9 already lead to the data compiled in Table |. 
We can note that three values lie below the prediction of Eq. (P , thereby confirming the semiclassical results given 
above and indicating that true degeneracy might not exist at £^ = for states of different Vs, after all. 

IV. CONCLUSION 

In conclusion, wc have shown that Eq. (|l|), proposed by Portnoi and Galbraith to describe the critical screening 
lengths for the two-dimensional statically screened Coulomb potential Vg, is unfortunately not strictly exact. This 
does not invalidate the variable-phase method approach, which is quite useful: the discrepancy between the exact 
values and the proposed integers is rather small. The reason for such a result is probably that Vg has the same 
asymptotic behavior than another potential, for which the critical lengths are exactly given by Eq. (|^) This 
should stimulate further work on two-dimensional cxcitonic screening, which is a subtle and exciting problem. 
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TABLE I. Upper bounds for (rs)^, for the lowest-lying states with I > 1. Values already below the prediction of Eq. (Q) are 
typed in bold face. 
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N=0 


N=l 


N=2 


N=3 


N=4 


N=5 


1 


3.003789 


5.98884 


10.0014 


15.1442 


21.6273 


29.6652 


2 


10.019610 


15.03181 


21.1347 


28.3099 


36.3503 


45.2024 


3 


21.034882 


28.00781 


36.0793 


45.3892 


56.1326 


69.0984 


4 


36.064718 


45.00193 


54.9819 


66.0182 


78.4080 


93.0400 


5 


55.107383 


66.02058 


77.9891 


91.6327 


105.8805 


122.6894 
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